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Abstract 

We calculate in baryon chiral perturbation theory the relativistic 1/M^-corrections to 
the leading order two-pion exchange diagrams. We give explicit expressions for the corre- 
sponding one-loop NN-amplitudes in momentum space. The resulting isovector central and 
isoscalar spin-spin and tensor NN-amplitudes involve non-static terms proportional to the 
. squared nucleon center-of-mass momentum p^. From the two-pion exchange box diagrams 

\ we obtain an isoscalar quadratic spin-orbit NN-amplitude. We give also analytical expres- 

CN ' sions for the corresponding NN-potentials in coordinate space. The diagrammatic results 

Q^. presented here make the chiral NN-potential complete at next-to-next-to-next-to-leading 

^ ! order. 

^ ■ PACS: 12.20.Ds, 12.38.Bx, 12.39.Fe, 13.75.Cs 
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Over the past years effective field theory methods have been successfully applied to the 
two- nucleon system at low and intermediate energies [|1], 0, ^, ||, ||, 0. The idea of constructing 
Q^ ' the NN-potential from effective field theory was put forward by Weinberg and this was first 
^ ■ taken up by van Kolck and collaborators [0 who used "old-fashioned" time-ordered perturbation 
Q . theory. Later, the systematic method of unitary transformations was employed by Epelbaum, 
^ ! Glockle and MeiBner |^ to construct an energy-independent NN-potential from the effective 
chiral pion-nucleon Lagrangian at (next-to-)next-to-leading order. Based on one- and two-pion 
CJ \ exchange and nine adjustable NN-contact interactions which contribute only to S- and P-waves 
^ ■ a good description of the deuteron properties as well as the NN phase-shifts and mixing angles 
below Tiab = 300 MeV was found in that framework 0]. Also recently, the elastic proton-proton 
scattering database below 350 MeV laboratory kinetic energy (consisting of 1951 data points) 
^ \ has been analyzed in terms of Ivr-exchange and chiral 27r-exchange at next-to-next-to-leading 
I order in ref.|0]. The resulting good x^/dof < 1 constitutes a convincing proof for the presence 
of the chiral 27r-exchange in the long-range proton-proton strong interaction. It was concluded 
in ref . that Ivr-exchange together with chiral 27r-exchange gives a very good NN-force at least 
as far inwards as r = 1.4 fm internucleon distance. All shorter range components of the NN- 
interaction have been effectively parametrized in ref.0 by 23 boundary condition parameters. 

At present, there is much interest in extending the calculations of the two- nucleon system 
[§] (and also the analysis of the NN database) to one order higher in the chiral expansion. 
This requires the full knowledge of the next-to-next-to-next-to-leading order (N^LO) chiral NN- 
potential. In momentum space this corresponds to terms in the NN T-matrix which are of 
fourth order in small external momenta and the pion mass, denoted generically by 0{Q'^). 
In particular, all two-loop diagrams of the process NN — > A^A^ with leading order vertices 
contribute at that order, 0{Q^). The chiral 3-7r-exchange has recently been calculated completely 
in ref. 0. Furthermore, results for the two-loop 2-7r-exchange are now available from ref. |T^. In 
that work also the one-loop contributions (of order 0{Q'^)) proportional to the second and third 
order low-energy constants Ci^2,3,4 and dj, related to elastic vrA-scattering, have been evaluated. 
The only remaining and presently unknown pieces of the NN T-matrix at order O(Q^) are the 
relativistic 1/M^-corrections to the leading order one-loop 27r-exchange. The purpose of this 
work is to present analytical results for these remaining contributions and thus to complete the 
chiral NN-potential at next-to-next-to-next-to-leading order. 



Let us first give some basic definitions in order to fix our notation. We are considering elastic 
NN-scattering in the center-of-mass frame, i.e. the process Ni{p) + N2{—p) — > Ni{p + q) + 
N2{—p — q). The corresponding on-shell T-matrix has the following general form, 



NN 



Vc + Ti- T2WC + [Vs + Ti ■ f2Ws\ (?i ■ ^2 + Wt + n ■ f2W^r] 0^1 • qa2 ■ q 

+ [Vso + Ti ■ ^"21^50] + ^2) ■ {qx p) 
+ fi • T2Wq] ai ■ (g X p) 0^2 ■ (g x p) , 



(1) 



with q = \q\ the momentum transfer between the initial and final-state nucleon, subject to the 
constraint 2p- q = —q^. The subscripts C, 5, T, SO and Q refer to the central, spin-spin, tensor, 
spin-orbit and quadratic spin-orbit components, each of which occurs in an isoscalar {yc,...,Q) 
and an isovector {Wc,...,q) version. Furthermore, ^1^2 and Ti^2 in eq.(l) denote the spin and 
isospin operators of the two nucleons. 



Fig.l: 2Tr-exchange bubble diagram. The combinatoric factor of this graph is 1/2. 



Let us now turn to the results for the 1 /M^-corrections to the chiral 27r-cxchange. The bubble 
or football diagram involving two isovectorial Tomozawa- Weinberg vertices is shown in Fig. 1. 
A convenient way to obtain all 1/Af ^-corrections from such a one-loop diagram goes as follows. 
At the left nucleon line one expands in the NN center-of-mass frame the pertinent relativistic 
TTTrNN-vertex sandwiched between Dirac-spinors for the out- and in-going nucleon in powers of 
1/M up to quadratic order, i.e. 0{M~'^). The analogous expression at the right nucleon line is 
obtained simply by switching the signs of three-momenta p and q and of the loop four-momentum 
{IqJ/}. Of course, one also relabels the spin and isospin operators, c?i (J2 and fi T2. In a 
last step, one multiplies the product of these both expressions with the pion propagators and 
carries out the four-dimensional loop integration. The total result for the 1/M^-corrections to 
the bubble diagram in Fig. 1 reads (displaying only non- vanishing contributions) , 

^ 7687r2M2/^ ' 

with the frequently occuring logarithmic loop function 

L(g) = ^ln^, ^ = ^4m^ + g2. (4) 
q zm^r 

In eqs.(2,3) and in the following wc omit piircly polynomial terms which can be absorbed in the 
strength of some zero-range NN-contact interactions. In fact, we are interested here only in non- 
polynomial or finite-range terms from chiral 27r-exchange. Note also that the isovector central 
NN-amplitude Wc in eq.(2) has a non-static contribution proportional the squared nucleon 
center-of-mass momentum p^ . 



Fig. 2: 2ti- exchange triangle diagrams. 



Fig. 2 shows the 27r-exchange triangle diagrams. The expression for the vrvrA^A^-vertex at the 
left (or right) nucleon line of the first (second) diagram can be taken over from the calculation of 
the bubble diagram. At the other (right or left) nucleon line one expands in powers of 1/M the 
product of upper pseudovector vrA^-vertex, nucleon Dirac-propagator and lower pseudovector 
vrA^-vertex sandwiched between out- and in-going Dirac-spinors. Collecting all possible terms 
proportional to and performing the loop integration over the pion propagators, the total 

result for the 1/M^-corrections to the triangle diagrams in Fig. 2 reads 



967r2M2/^ 



|— g"^ + 5m^g^ + 3m^ — 6m^w ^ — p'^^Sml + 5g^)| , (5) 



glLiq) , 3 



Again, the isovector central NN-amplitude Wc in eq.(5) has a non-static piece proportional to 
p^. In the heavy nucleon mass limit M — oo the diagrams in Fig. 1 and Fig. 2 give rise only to 
a non- vanishing isovector central NN-amplitude Wc (see eq.(14) in ref.|^) and one finds that 
the non-static p^-terms in eqs.(2,5) are in fact proportional to it. 




Fig. 3: 2t[ -exchange box diagrams. The (right) planar box graph includes the iterated -exchange. 



Fig. 3 shows the 27r-exchange crossed box and planar box diagrams. The crossed box diagram 
(of isospin structure 3 + 2ri ■ T2) can be evaluated straightforwardly by making use of the 
expression for the sequential pion absorption and emission process derived in the calculation of 
the triangle diagrams. At the left and the right nucleon line these expressions differ just by the 
sign of the three-momenta p and q. The planar box diagram (of isospin structure 3 — 2fi ■ 
requires some special treatment since it contains the so-called iterated Ivr-exchange. As outlined 
in section 4.2 of ref. one must perform the dl^ loop integration (using e.g. residue calculus) 
before the 1/M-expansion in order to isolate the iterated Ivr-exchange together with its own 
relativistic 1/M-corrections. From that (and also from a consideration of the Lorentz- invariant 



NN two-body phase space) one finds a dependence of tlie iterated Ivr-exchange on tlie nucleon 
mass M of the form M'^/y/M'^ = M - p^/2M + O^M'^). Since there obviously is no 
contribution to the iterated Ivr-exchange proportional to one can treat (for the present 

purpose) the planar box graph in the same way as the crossed box graph. Namely, one performs 
first the 1/M-expansion of vertices and propagators at the left and the right nucleon line and 
then one carries out the four- dimensional loop integration. We have used this method as well 
as the continuation of the method outlined in section 4.2 of ref.0 (where things are basically 
done in the reverse order) in order to calculate the 1/M^-corrections to the planar box graph 
and we have found perfect agreement. Sorting finally the contributions from the crossed box 
diagram and planar box diagram in Fig. 3 into isoscalar and isovector components, we find after 
a tedious calculation 

Vc = 32jlr2.4 {^<^"' + [2<^"' + 8<^"' - - 2<]L(g)} , (8) 
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Wc = — Hm^'w 



27 



L{q)}, (9) 



1... 9\H(l) (^2 , 17 2 , .4. .-2 



9aH(1) fll„2 , _4„,-2 



Vn = -^^^. (14) 

Here, one observes that the isovector central NN-amplitude Wc and the isoscalar spin-spin 
and tensor NN-amplitudes Vs^t in eqs.(9,10) receive non-static contributions proportional to p^. 
Another interesting feature is the occurrence of an isoscalar quadratic spin-orbit NN-amplitude 
Vq in eq.(14), whereas the analogous isovector NN-amplitude vanishes, Wq = 0. This completes 
the presentation of analytical results for the momentum space NN-amplitudes. 

For certain applications 0] it is useful to have a representation of the NN-interaction in 
terms of local coordinate space potentials. These potentials are related to the momentum space 
NN-amplitudes by a negative Fourier-transform, — (27r)"^ / d^q exp{iq-f). The isoscalar central, 
spin-spin, tensor and spin-orbit potentials, denoted here by Vc,s,T,so{r), are accompanied by 
the operators 1, ?i ■ ^2, 3 o^i ■ f c?2 ■ — ?i ■ ?2 and — |(c?i + c?2) ■ (r x V). In the case of the 
isovector potentials Wc,s,T,so{r) an additional isospin operator ri ■ T2 occurs. The finite-range 
parts (disregarding zero-range (5^(r )-terms) of the coordinate space potentials following from 
eqs.(2-14) are most conveniently obtained through their spectral-representations p|, |T0| with the 
help of the formula Im L{ifi) = —{7i/2^)^n'^ — Am'^. It turns out that all these potentials can be 
written analytically in terms of two modified Bessel-functions Kq i{2x), with the dimensionless 
variable x = m.,^r. We list them separately for the three classes of diagrams. 

i) Bubble diagram in Fig. 1: 

Wc{r) = i287r3M2/4r6 {^"'^^ ^ ^ p'r^)M2x) + (10 + 7x^ + 2p'r^)K,{2x)} , (15) 



^'^'^ ^ 3847r3MV.v {^'^^^ + a;2)A'o(2a;) + (10 + lx')K,{2x)] , (16) 

^^^'^ ^ - 1536 JmV.V^ ^ ^x^)K,{2x) + (35 + 20x^)K,(2x)} , (17) 

^^^^"^ ^ "256;S^{^'^^°^^^^ + + 2x^)X,(2x)} . (18) 

ii) Triangle diagrams in Fig. 2: 

^c{r) = ^^J^HllO + 34a;^ + a:^+pV)iro(2a;) 

+[110 + 89x^ + + (10 + 4x^)pV]Xi(2x)} , (19) 

^'^'^ ^ - 192^^^^,6 (10^(7 + 2x^)i^o(2x) + (70 + 55x^ + Ax')K,{2x)} , (20) 
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^^^"■^ = 768^31^^2/^6 {^(245 + 52x^)Xo(2x) + (245 + 170x^ + 8x')K^{2x)} , (21) 

^so{r) = i28JwV>^ {^^^^ + 4x^)Xo(2x) + (45 + 2ix')K^{2x)} . (22) 

iii) Box diagrams in Fig .3: 

^cir) = - isJmV^^ I^^^" + ^^"^^ + x'Wx) + (30 + 27a;^ + ^x' + ^)i^i(2a;)} , (23) 



^^^""^ ^ agJaWre t"^^'^'^^"^' - 42) - 810 - 258x-^ - 6x"]Ko{2x) 

-[810 + 663x^ + 70j;^ + 4x^ + ^^^(42 + Sx'^)]Ki{2x)} , (24) 



^'^'^ ^ T^^^?^H90 + 36^' + 4^'-12pV)^o(2x) 

+[90 + 81^2 + 14x^ - pV(12 + 8x'^)]Ki(2x)} , (25) 



^5(r) = /^Tl; , 6 (^(510 + 162x2 + 8x^)iro(2x) + (510 + 417^^ + Ux^)Ki{2x)} , (26) 
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^^^'^ ^ 768'3MV4,a {^(48A^ - - HAx^ - 8x^)i^o(2x) 



+ [pV2(60 + 16^2) - 315 - 270x2 - 40a;^]is:i(2a;)} , (27) 

4 

VrT(r) = - 46og^1^2j4,6 {^(17S5 + 456^2 + 16x*)Ko{2x) + (1785 + 1320^^ + n2x*)K,{2x)} , 

(28) 

^so{r) = IsS^^TM^"^^^^^ + 52x2)i^o(2x) + (165 + 132^2 + 16x^)i^i(2x)} , (29) 
^^^^^^ ^ 768^W^ " "'^^^^ ^ 4x2)Ko(2x) + (8x^ - 78^^ - 165)K,{2x)} , (30) 



= ,.it7f4.6 {l2a:i^o(2x) + (15 + Ax')K,{2x)] , (31) 



647r3M2/4r6 

The last potential Vqlr) is accompanied by the quadratic spin-orbit operator —ai ■ (r x V) 0^2 • 
(r X V). In order to Fourier-transform the w~^L(g)-terms in eqs.(8,9) it is useful to exploit the 
relation 4m'^w~^L{q) + w'"^ = — m^(9/9m^)[w~2L(g)]. 

In Table 1, we present numerical values for the coordinate space potentials eqs. (15-31) at 
p=0, using the parameters gA = 1-3, = 138 MeV, U =^2.4 MeV and M = 939 MeV. One 
observes that the attractive central potentials, Vc{r) and Wc{r), and the repulsive isoscalar 
spin-orbit potential Vso{r) are quite sizeable. The inclusion of the relativistic l/M^-corrections 
to the chiral 27r-exchange in NN phase-shift calculations etc. may therefore not be just a small 
effect. Of course, a firm conclusion about their relevance can only be drawn from the complete 
N^LO chiral NN-potential. 

In summary, we have completed here the chiral NN-potential at next-to-next-to-next-to- 
leading order by evaluating the relativistic l/M^-corrections to 27r-exchange. The analytical 
results presented here are in a form such that they can be easily implemented in a N^LO 
calculation of the two-nucleon system or in an empirical analysis of low-energy elastic NN- 
scattering. Work along this line is in progress |p. With such calculations one will be able to 
determine the maximal range of validity of chiral multi-pion exchange in the NN-interaction. 



r [fm] 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


Vc [MeV] 


-185.5 


-72.3 


-31.4 


-14.9 


-7.56 


-4.06 


-2.29 


-1.34 


Wc [MeV] 


-103.7 


-40.2 


-17.4 


-8.18 


-4.13 


-2.20 


-0.28 


-0.72 


Vs [MeV] 


46.4 


18.1 


7.86 


3.73 


1.89 


1.02 


0.57 


0.34 


Ws [MeV] 


23.1 


8.95 


3.87 


1.83 


0.92 


0.49 


0.28 


0.16 


Vt [MeV] 


-40.2 


-15.6 


-6.77 


-3.20 


-1.62 


-0.87 


-0.49 


-0.29 


Wt [MeV] 


-19.9 


-7.69 


-3.32 


-1.56 


-0.78 


-0.42 


-0.23 


-0.13 


Vso [MeV] 


124.8 


48.3 


20.9 


9.81 


4.94 


2.63 


1.47 


0.85 


Wso [MeV] 


-2.15 


-0.78 


-0.31 


-0.13 


-0.06 


-0.028 


-0.013 


-0.006 


Vq [MeV] 


19.4 


7.26 


3.02 


1.37 


0.66 


0.34 


0.18 


0.10 



Tab.l: Numerical values of the relativistic 1/ M"^- corrections to the chiral 2'K-exchange NN- 
potentials (at p = 0) versus the nucleon distance r. The units of these potentials are MeV. 
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